Observational Constraints on finite scale factor singularities 



O 
(N 

^3 



Tomasz Denkiewicj^ 

Institute of Physics, University of Szczecm, Wielkopolska 15, 70-451 Szczecin, Poland and 
Copernicus Center for Interdisciplinary Studies, Slawkowska 17, 31-016 Krakow, Poland 

(Dated: June 28, 2012) 

We discuss the combined constraints on a Finite Scale Factor Singularity (FSF) universe evolution 
scenario, which come from the shift parameter TZ, baryon acoustic oscillations (BAO) A, and from 
the type la supernovae. We show that observations allow existence of such singularities in the 2 x 10^ 
years in future (at la CL) which is much farther than a Sudden Future Singularity (SFS), and that at 
the present moment of the cosmic evolution, one cannot differentiate between cosmological scenario 
which allow finite scale factor singularities and the standard ACDM dark energy models. We 
also show that there is an allowed value of m = 2/3 within la CL, which corresponds to a dust- 
filled Einstein-de-Sitter universe limit of the early time evolution and so it is pasted ito a standard 
early-time scenario. 
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I. INTRODUCTION 

Finite scale factor singularities (FSF) are one of types 
of new and exotic singularities which were found first in 
Ref. [l|. The inspiration to search these new types of 
singularities was due to the observations of high-redshift 
type la supernovae (SNIa) which provided strong evi- 
dence that the expansion of the universe is accelerating 
due to an unknown form of dark energy [s^ phenomeno- 
logically behaving as the cosmological constant. Further 
observational data [1] made cosmologists think of more 
accelerating universe filled with phantom ^ which vi- 
olated all energy conditions: the null [qc^ -f p > 0), 
weak {qc^ > and gc^ + p > 0), strong {gc^ + p > 
and gc^ -\- 3p > 0), and dominant energy {gc^ > 0, 
—gc^ < p < gc^) (c is the speed of light, g is the mass 
density in kgm~^ and p is the pressure). Phantom-driven 
dark energy leads to a big-rip singularity (BR or type I 
according to Q ) in which the infinite values of the energy 
density and pressure (p, p — oo) are acconipanied by the 
infinite value of the scale factor [a — > oo) |6|. 

The list of new types of singularities contains: a big- 
rip (BR), a sudden future singularity (SFS) which 
can appear in inhomogeneous and anisotropic models too 
a generalized sudden future singularity (GSFS), a 
finite scale factor sing ularity (FSF) 'IS^, a big-separation 
singularity (BS). They are characterized by violation of 
all or some of the energy conditions which results in a 
blow-up of all or some of the physical quantities: the scale 
factor, the energy density and the pressure. The finite 
scale factor singularity, which is the subject of this pa- 
per, is a weak singularity according to Tipler's definition, 
but strong according to Krolak's definition [13]. Apart 
from mentioned above there are also w-singularities |15| , 
which are not physical singularities, but are singularities 
of a barotropic index w which are present in different 



cosmological models of f{R) gravity [16], in scalar field 
models [iTj, and in brane cosmologies [18] . 

For a finite scale factor singularity (FSF) p — >■ oo and 
p — > oo diverge, while the scale factor a remains constant. 
This means that it is similar to a big-bang singularity 
with only one exception - the scale factor is a constant 
a = as = const instead of zero ass = 0. Besides, FSF 
singularities are stronger than a sudden future singular- 
ities (SFS) [10], so they placed themselves in between 
the big-bang (which is strong and geodesically incom- 
plete) and the SFS. It has been shown that SFS appear 
in physical theories such as f{R) gravity, scalar field cos- 
mologies, brane cosmologies [191 ] and, in particular, they 
plague loop quantum cosmology f23|. It seems then that 
after an appropriate choice of the scalar field, f{R) func- 
tion or brane, FSF will be easily deduced to appear in 
such physical theories as well. 

The paper is organized as follows. In section [ll] we 
present FSF scenario. In section lllll we derive the expres- 
sions for the observables: baryon acoustic oscillations, 
distance to the last scattering surface, type la super- 
novae, used to test an FSF scenario. In section IIVI we 
give the results and discussion. 



II. FINITE SCALE FACTOR SINGULARITY 
SCENARIO CONSTRUCTION 

In order to obtain an FSF singularity one should con- 
sider the simple framework of an Einstein- Friedmann cos- 
mology governed by the standard field equations 
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is trivially fulfilled due to the Bianchi identity. Here a = 
a{t) is the scale factor, the dot means the derivative with 
respect to time i, G is the gravitational constant, and the 
curvature index fc = 0, ±1. During further considerations 
we set c — \ and fc = 0. 

Similarly like in the case of SFS, which were tested 
against the observations in [2l| - |23| . one is able to obtain 
an FSF singularity by taking the scale factor in the form 

a[y) = a,[5+{l-5)jr-S{l~yT] , y^f 

(11.4) 

with the appropriate choice of the constants (5, ts,as, m, n. 
In contrast to an SFS in order to have accelerated ex- 
pansion of the universe 6 has to be positive (5 > 0). For 
1 < n < 2 we have an SFS, which plagues the loop quan- 
tum cosmology [20| . In order to have an FSF singularity 
instead of SFS, n has to lie in the range < n < 1. 

As can be seen from (jII.l|) - (|II.4p . for an FSF p diverges 
and we have: for t ^ ts, a ^ a^, p — > oo, and \p\ — > oo, 
where Os, tg, Ps, are constants and Us ^ 0. 

In a model, expressed in terms of the scale factor (|IL4p , 
the evolution begins with the standard big-bang singular- 
ity at t = for a = 0, and finishes at a finite scale factor 
singularity for t = tg where a = = a{ts) is a constant. 
In terms of the rescaled time y we have a(l) = as. 

The standard Friedmann limit (i.e. models without a 
singularity) of (jll.4l) is achieved when (5 — 0; hence 6 is 
called the "non-standardicity" parameter. Additionally, 
notwithstanding Rcf. [7] and in agreement with the field 
equations pLT|) - (llL2]) . S can be both positive and neg- 
ative leading to an acceleration or a deceleration of the 
universe, respectively. 

It is important to our discussion that the asymptotic 
behaviour of the scale factor (jll.4[) close to the big- 
bang singularity at t = is given by a simple power- 
law flBB = 2/"*, simulating the behaviour of flat k = 
barotropic fiuid models with m — 2/[3{'w + 1)] . 

The FSF singularity scenario consists of two dominat- 
ing components such as a nonrelativistic matter, and the 
fiuid which is driving the singularity. We consider the 
case of the noninteracting components from which both 
of them obey independently their continuity equations as 
III. 31 The evolution of both ingredients is independent. 
Nonrelativistic matter scales as 

—j (11.5) 

and the evolution of the other fluid, which we name here 
Pq, can be determined by taking the difference between 
whole energy density, p, evolution from Friedmann eq. 
III. II and the pm- 

PQ= P- Pm (11.6) 

this ingredient, pg, of the content of the Universe is re- 
sponsible for the singularity in the energy density for 

t ts. 



A. Classical stability of FSFS with respect to small 
perturbations 

In Ref. the study of classical stability of sudden 
singularities and "Big Rip" singularities with respect to 
small inhomogeneous scalar, vector and tensor perturba- 
tions was given. The analysis was made in the framework 
of g auge-invariant formalism introduced by Mukhanov 
|10| . Characterization of sudden singularities in terms of 
the series expansion of the scale factor on the approach to 
a singularity was given, and it was shown that once the 
density diverges near to a singularity it is unstable under 
small scalar perturbations. On the other hand, sudden 
singularities are stable to inhomogeneous vector and (ten- 
sor) gravitational-wave perturbations. Taking the Fried- 
mann equations (|II.ip - (jII.2[) with c — 1, 8ttG — 1 and 
the scale factor in terms of the series expansion on the 
approach to a singularity: 

a{t)=co + ci{U-t)^ + ..., (II.7) 



a{t) ^ -ciX{t, - ty-^ + . . . , (II.8) 



d{t) = ciA(A - l){ts - tf-^ + ..., (II.9) 

were the c^, A are real constants, with Ci 7^ 0, cq > 0, in 
the limit t ^ ts the approximation for the density and 
pressure was derived [9|: 

p=^ + ^(t,-t)2(A-l)+..., (11.10) 



p^-^^'-^^i^^it.^tr^^... (11.11) 

Cq Co 

The density diverges if and only if A < 1, and the 
pressure diverges, if A < 2. Substituting pi.7IIII.lll) into 
the gauge- invariant scalar perturbations equation [10|: 

a'^<i> + {4: + 3p{p))aa^ + 
{2a& + (a^ - fc)(l-f^ 3p'{p)) + p'{p)K^)<^> = 0, (11.12) 

for plane wave perturbations with wave number K, where 
p'{p) = dp/ dp, and further neglecting higher-order terms, 
and taking the case of A < 1 we have j9|]: 

$ - (A - 2)r"i<i> + — A^r^"^^ = 0, (11.13) 
Co 

where T = — t and dots indicate differentiation with 
respect to t. By transformation of the variables, one can 
express this equation as the Bessel equation, and then 
obtain its solution to be H: 
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B. Shift parameter 



(11.14) 



where A, B are constants. For FSFS with A < 1, 
diverges as t — >■ ts- Such a result is vahd for all 
wavelengths, and for all values of k. 



III. OBSERVATIONAL CONSTRAINTS 

We consider three observational constraints, such as lu- 
minosity distance moduli to type la supernovae, baryon 
acoustic oscillations, and the shift parameter which is 
scaled distance to the last scattering surface of the cos- 
mic microwave background. We make use of Markov 
Chain Monte Carlo technique within the framework of 
Bayesian statistics to obtain posterior probability distri- 
bution as a function of the parameters: n, m, (5, j/q = 
where is the present age of the universe. We em- 
ploy Metropolis-Hastings algorithm with uniform priors: 
5 e (0, 30), n e (0, 1), TO e (0, 3), e (0, l). The results 
are marginalised over Hq = 74.2 ± 3.6 taken from HST 
Key Project. 



A. Supernovae 

We proceed within the framework of Friedmann cos- 
mology, and consider an observer located at r = at 
coordinate time t = Iq. The observer receives a light ray 
emitted at r = ri at coordinate time t = ti. We then 
have a standard null geodesic equation 
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(III.l) 



with the scale factor a{y) given by (|II.4|) . Using (jll.4p 
again, the redshift is given by 

a{h) 6+{l-S)yr -S{l-yir ^ ' 

where yo = y{to) and j/i = y{ti). 

The luminosity distance to supernova is given by 



dL{z) = (1 + z)cts 



dy 



and the distance modulus is 

/i(z) = Slogio dL{z) + 25. 
The for the SNIa data is 

2 _ \^ (Aio6s(2i) - 
XSN 



1=1 



(III.3) 



(III.4) 



(in.5) 



where ai is the quoted observational error on the z**^ 
Union2 SNIa and CTint is the SNIa intrinsic scatter. Fol- 
lowing |2J| we take crint =0.15. 



The standard formula for the CMB shift parameter is 
given by [H ill: 



7^ = 



7/TT 

n 

iTT 



(III.6) 



where ij"^ is the temperature perturbation CMB spec- 
trum multipole of the first acoustic peak in the model 
under consideration and I'l"^ corresponds to a reference 
fiat standard Cold Dark Matter (CDM) model. The mul- 
tipole number is related to an angular scale of the sound 
horizon at decoupling by [13, [llj 



dA h 



where 



Ts = adccS{rs) = adccS ( / '^^^^ 
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with Cs being the sound velocity and the angular diameter 
distance reads as 



= adoc<5'(?'dcc) = adcc'S' ( / c-^ 

tdcc «(^) 
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where c is the velocity of light and S{r) — r for fc = 0. 

Following [2^ we can write, using (|III.6p and (|III.7|) as 
follows 



n = 
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which, by assuming that at decoupling the amount of 
radiation was the same in both the flat reference standard 
CDM model and in our FSF model (which we assume to 
be just the same as a standard matter-radiation model of 
the early universe, since FSF models do not change the 
evolution there) we have that 
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On the other hand, for a reference standard CDM model 
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while for our FSF model the angular diameter distance 
is given by 



dA 



(III.13) 



with r^cc given by taken at decoupling. Using the above, 
we may write that for our FSF models the shift parameter 
is 

2c 2ct, 



n 



Haaoy/Qmordcc a'{yo)\/^mQrdc 



(III.14) 
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where we have assumed that the function /(fl^oi '^dcc) 
approximately unity [2^ . 

Finally, the rescaled shift parameter is 
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where E{z) = H{z)/Hq. The WMAP data gives 7^ = 
1.725 ± 0.018 [H], thus ^oi the shift parameter is the 
following 



2 _ (7^- 1.725)2 
Xr - 



0.0182 



C. Baryon acoustic oscillations 



(III.16) 



In Ref. [30j it was suggested that instead of taking 
the position of an acoustic peak one should measure the 
large-scale correlation function at lOOh^^ Mpc separation 
using the sample of 46748 luminous red galaxies (LRG) 
selected from the SDSS (Sloan Digital Sky Survey) sam- 
ple. The appropriate quantity to be measured is known 
as the A parameter and reads as 



A = nli^EizBAoY 
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In Fig. [T]we see that there are two disjoined regions 
in the parameter space for which FSF singularities are 
allowed in the universe. 6 is always positive which is 
not surprising since only for positive S there can be an 
accelerated expansion in the investigated model. Char- 
acteristic feature of the picture is that there are two qual- 
itatively different regions divided by the value of the pa- 
rameter m. There are two branches, the first for m > 1, 
and the second for m < 1. In the following tables we 
present ranges of the values of the parameters for three 
confidence levels. For the case when m > 1: 



la CL 


2a CL 


3a CL 


m e (1.00, 1.14) 


m e (1.00,1.21) 


TO e (1.00,1.28) 


5 > 1.00 


S > 1.00 


S > 1.00 


n e (0.0,0.97) 


n e (0.0,0.98) 


n e (0.0,0.99) 


yo e (0.0,0.50) 2/0 e (0.00,0.66) yo G (0.00,0.76) 


Below the case of 


TO < 1: 




la CL 


2a CL 


3ct CL 


m £ (0.41,0.66) 


m e (0.39,0.67) 


TO e (0.39,0.72) 


6 € (0.52,0.99) 


6 e (0.36,0.99) 


5 G (0.26,0.99) 


n e (0.0,0.85) 


n e (0.0,0.91) 


n e (0.0,0.94) 


2/0 e (0.45,0.88) 


yoe (0.43,0.94) yo e (0.41, 0.97) 



According to Ref. [30| A should have the value 



A = 0.469 



Vo.98/ 



-0.35 



±0.017 



(III.18) 



where n is the spectral index (now taken about ^ 0.96). 
We have for BAO 



2 (yl- 0.469)2 
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0.0172 



(III.19) 



Using pL4)) and (|nL2t we get the form of A for an FSF 
as follows 
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IV. RESULTS AND CONCLUSIONS 



Overall that has been used is 

X^=xIn+xI + Xa (IV.l) 

In Fig. [T] we present contour plots showing the joint 
marginal posterior distribution for each pair of FSF 
model parameters. Each sub-panel shows three contours, 
denoting roughly 68%, 95% and 99% (from light gray to 
dark grey respectively) credible regions. 
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FIG. 1: Marginalized contours for pairs of parameters are plotted. There are three confidence regions 68%, 95%, 99% (from 
hght grey to dark grey respectively) calculated for A, TZ, and SN la jointly 



For the case of the branch m > 1 the maximum value 
at 3(7 uq parameter is 0.76, which means that such a 
singularity can happen later than for the second branch 
m < 1. While yo approaches 0.76, S grows stronger and 
as the value of 6 more grows, m becomes better con- 
strained. 

What seems to be important is that there is an allowed 
value of rn = 2/3 within la CL, which could correspond 
to the dust filled Einstein-de-Sitter universe in a close to 
big-bang limiting case. What is also interesting for this 
branch is that while the allowed values of the nonstan- 
darcity parameter S are small, i.e. 6 < 1, the parameter 
Uo approaches unity which means that the singularity 
may happen in the nearest future for this case. For la 
CL, uq = 0.885 corresponds to ^ 2 x 10^ years to the time 
of the singularity. For 3a CL, yo = 0.97 corresponds to 
the present time at ~ 0.37 x 10^ years before the time of 
the singularity. 

In conclusion, we have shown that for a finite scale 
factor singularity there is an allowed value of m = 
2/3 within la CL, which corresponds to the dust-filled 



Einstein-de-Sitter universe for the close to big-bang lim- 
iting case. The finite scale factor singularity may happen 
within 2 x 10^ years in future for la CL and its predic- 
tion at the present moment of cosmic evolution cannot be 
distinguished, with current observational data, from the 
prediction given by the standard quintessence scenario of 
future evolution in the Concordance Model [3ll - l40| . 
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